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We discuss two models for describing the behavior of matter at large densities and intermediate 
temperatures. In both models a softening of the equation of state takes place due to the appearance 
of new degrees of freedom. The first is an hadronic model in which the softening is due to chiral 
symmetry restoration. In the second model the softening is associated with the formation of clusters 
of quarks in the mixed phase. We show that both models allow a significant softening but, in the 
first case the bulk modulus is mainly dependent on the density, while in the mixed phase model it 
also strongly depends on the temperature. We also show that the bulk modulus is not vanishing 
in the mixed phase due to the presence of two conserved charges, the baryon and the isospin one. 
Only in a small region of densities and temperatures the incompressibility becomes extremely small. 
Finally we compare our results with recent analysis of heavy ion collisions at intermediate energies. 

PACS numbers: 21.65.+f, 12.39.Fc, 25.75.-q, 26.60.+C. 
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Introduction. - The behavior of matter at large densi- 
ties and temperatures is still poorly known but, on gen- 
eral grounds, new degrees of freedom are expected to 
appear and they should generate a softening of the Equa- 
tion of State (EOS). In this Letter we discuss two mod- 
els for the softening: the first model is based on a purely 
hadronic chirally symmetric EOS and the soften- 

ing is due to partial restoration of the chiral symmetry. In 
that model the calculations are performed beyond mean 
field approximation and quantum fluctuations play a cru- 
cial role. In a second model the softening is due to the 
formation of clusters of quarks, which are the precursors 
of deconfinement. These clusters are at first metastable 
and they stabilize only at larger densities. 

The aim of our calculation is to compare how the soft- 
ening takes place in the two models and finally to relate 
our results to recent analysis of the experimental data. 
This may be helpful also in prevision of future experi- 
ments planned e.g. at facility FAIR at GSI 

The extraction of experimental information about the 
EOS of matter at large densities and temperatures from 
the data of intermediate and high energy Heavy Ion Col- 
lisions (HICs) is very complicated. Possible indirect indi- 
cations of a softening of the EOS at the energies reached 
at AGS have been discussed several times in the litera- 
ture d, 0, 0, H, H 113 . In particular, a recent analysis 
[m based on a 3-fluid dynamics simulation suggests a 
progressive softening of the EOS tested through HICs at 
energies ranging from 2A GeV up to 8A GeV. 

In this letter we will show that in the chiral model the 
incompressibility strongly depends on the density but it 
is almost temperature independent at least up to T 150 
MeV. At the contrary a strong temperature dependence 
is found in the mixed phase model. Since in HICs at in- 
termediate energies not too large densities arc reached, 
the mixed-phase model seems to provide a better descrip- 
tion of the results of the experiments analysis. 



The Chiral- dilaton model. - It is not trivial to develop 
chiral invariant models in which saturation properties of 
nuclei are well reproduced. It is well known that the 
attempt of using e.g. the sigma model to describe nuclear 
dynamics fails due to t^he impossibility to reproduce basic 
properties of nuclei [12]. More sophisticated approaches 
have been proposed in the literature, both within a SU(2) 
chiral symmetric model [l^, [l3| and also extending the 
symmetry to the strange sector [H, [l^, [13] ■ Here we 
use the model introduced by the Minnesota group [J, 
13, [3- that model chiral fields are present together 
with a dilaton field which reproduces at a mean field 
level the breaking of scale symmetry which takes place 
in QCD. In [l|, it has been developed a formalism 

(which we adopt) allowing resummations beyond mean 
field approximation. This is important when studying a 
strongly non-perturbative problem as the restoration of 
chiral symmetry. The lagrangian of the model reads: 
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Here a and tt are the chiral fields, the dilaton field, 
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FIG. 1: The left panel represents the free energy of a bubble 
of quarks as a function of its radius R, while in the right panel 
the temperature vs. the baryon density is shown. The a and 
e curves correspond to the critical densities delimiting the 
mixed phase in the absence of surface tension; the densities 
between b and c correspond to the formation of a metastable 
bubble. For densities larger than c the mixed phase is stable 
(an example of the corresponding free energy is provided by 
the d curve). 



the vector meson field and the vector-isovector me- 
son field, introduced in order to study asymmetric nu- 
clear matter. The field strength tensors are defined in 
the usual way F^^, = d^^Uy - dyOJ^, B^^ = d^h^ - 9^b^. 
In the vacuum (f) — (/jq, cr = ctq and tt = 0. The lu and 
p vacuum masses are generated by their couplings witli 

1/2 1/2 

the dilaton field so that TOo; = GJ^ 4>o and rrip = G^^ cj)o. 
Moreover C, — (po/ao, B and S are constants and e[ is a 
term that breaks explicitly the chiral invariance of the 
lagrangian. The potential V of Eq. ([2]) is responsible for 
the scale symmetry breaking. The choice of such a po- 
tential comes from the necessity to reproduce the same 
divergence of the scale current as in QCD. In our calcula- 
tion we use the parameters set of Ref. [3] which was able 
to reproduce nuclear spectroscopy and also gives, within 
this model, the smallest value for the incompressibility 
at saturation density, K^^ — 322 MeV. Notice that this 
value is slightly larger than those traditionally used. 

The mixed hadron-quark phase. - If the deconfinement 
transition at finite density is first order, a mixed phase 
can form and it is typically described using two separate 
EOSs, one for the hadronic and one for the quark phase. 
Concerning the hadronic phase we use a relativistic field 
theoretical model, the NLpi5 [l8[, taking into account also 
a scalar-isovector interaction, which increases the sym- 
metry energy only at large densities. Qualitatively sim- 
ilar results can be obtained using other hadronic mod- 
els. For the quark phase we adopt an MIT bag like 
model. It is well known that, using the simplest version 
of the MIT bag model, if the bag pressure B is fixed to 
reproduce the critical temperature computed in lattice 
QCD, then, at moderate temperatures, the deconfine- 
ment transition takes place at very large densities. On 
the other hand there are strong theoretical indications 
that at moderate and large densities (and not too large 
temperatures) diquark condensates can form, whose ef- 
fect can be approximately taken into account by reducing 
the value of the bag constant. A phenomenological ap- 



proach can therefore be based on a density dependent bag 
constant, as proposed in Refs.fTol. [20|. We have adopted 
a parametrization of the form 

Beff = Bo- [A(m)]'m' , (3) 

where A(^) = Aexp[-(/i - po'f/a'^]. Here Bq^^'^ = 215 
MeV, A = 100 MeV, po = 300 MeV and a = 300 MeV. 
One gluon exchange corrections are taken into account 
and we use as — 0.35. One constraint on the parameter 
values is that a,t p — the critical temperature is ~ 170 
MeV, as suggested by lattice calculations, while the other 
constraint is the requirement that the mixed phase starts 
forming at a density slightly exceeding 3 po for a temper- 
ature of the order of 90 MeV (as also suggested e.g. by 
[2l|). The choice of the parametrization is clearly in- 
spired by the results of microscopical analysis on diquark 
condensate (for recent reviews see e.g. [12, HI]). On the 
other hand, the actual formation of a condensate and, 
even more, its specific type are still very uncertain [ssj . 
what justifies a phenomenological approach as ours. 

To describe the mixed phase we use the Gibbs formal- 
ism, which in Refs. [l^, [IE [ll] has been applied to sys- 
tems where more than one conserved charge is present. 
In this Letter we are studying the formation of a mixed 
phase in which both baryon and isospin charge are pre- 
served [36| . The main result of this formalism is that, 
at variance with the so called Maxwell construction, the 
pressure is not constant in the mixed phase and therefore 
the incompressibility does not vanish [STf . 

An important issue concerns the effect of a surface ten- 
sion at the interface between hadrons and quarks. The 
value of such a tension is poorly known, and in an MIT- 
bag-like model it is dominated by the effect of finite 
masses J2^, in particular of the strange quark. On the 
other hand in the system studied here strangeness plays 
a minor role, if any, and therefore the surface tension a 
should be rather small. In the analysis we have used a 
— 10 MeV/fm^, but the results are qualitatively similar 
if a slightly larger value of a is used. In the left panel 
of Fig. [1] we show the free energy for the nucleation of a 
bubble of quarks for various values of the density, com- 
puted using the formalism of Ref. [1^. One can notice 
that the density region in which quark bubble nucleation 
can really form shrinks with respect to the case with 
(7 = 0. Moreover, there is an intermediate density re- 
gion in which bubbles can form, but they are metastable. 
Only at larger densities a stable mixed phase can form. 
In our calculation this region is rather tiny, as shown in 
the right panel of Fig. [U but its size would increase for 
larger values of a. This region of metastable quark bub- 
bles is conceptually very interesting, because it provides 
a link with the so-called hadron-string-dynamics model 
[lH- In both descriptions, metastable matter is produced 
in an intermediate density window. In the string forma- 
tion scenario that matter is made of unstable hadrons 
while in the scenario here presented it is made of (small) 
unstable bubbles of quarks. 
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FIG. 2: Pressure at T=0 for symmetric nuclear matter, ob- 
tained using the chiral-dilaton model, the NLp5 model, the 
GM3 model and the set of simple parametric EOSs used in 
Ref. '11] . Also the pressure obtained using the GM3 with a 
small hyperon fraction is shown. 
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Results - We start by comparing in Fig. [5] the pres- 
sure computed using the models here discussed with the 
limits obtained from the analysis of HICs at intermedi- 
ate energies [30]. We have also introduced an EOS with a 
fixed fraction of strangeness due to the presence of hyper- 
ons (using the GM3 parametrization). In HICs a small 
number of hyperons is generated through associated pro- 
duction. An estimate of the number of hyperons per 
participant can be obtained from the experimental ratio 
of the kaon yields per participant sT], which indicates 
a strangeness fraction smaller than 10% _38]. As a ref- 
erence, we also show the parametric EOSs discussed in 
[ll| . which are characterized by different values of the 
incompressibility parameter K^^. From Fig. [2] it is clear 
that only the extra soft and the hard EOSs are excluded 
by the experimental limit. The pressure computed us- 
ing the chiral-dilaton model marginally exceeds the limit 
at low densities, due to the too large value of the in- 
compressibility at saturation density. The effect of the 
partial restoration of chiral symmetry is clearly visible as 
a softening taking place at larger densities. Instead the 
small fraction of strangeness introduced using the GM3 
parametrization does not produce a sizable softening of 
the EOS. 



FIG. 3: The upper and lower panels show respectively the pa- 
rameter and the bulk modulus at different temperatures 
as a function of density for the chiral-dilaton model and the 
GM3 parametrization. Here we used Z/A = 0.4. 
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FIG. 4: Same as in Fig. 3, but using the mixed-phase model. 



We now compare the compressibility computed in the 
various models discussed above with the one estimated in 
[ll| . In the following we always assume that Z/A = 0.4. 
In the lower panels of Figs. [3] and |4] we plot the bulk 
modulus B = pdP/dp, respectively for the chiral dila- 
ton model and for the mixed-phase model. Similarly, in 
the upper panels we show, as a function of the density 
and temperature, the value of the incompressibility pa- 
rameter for which the bulk modulus of a parametric EOS 
has the same value of the bulk modulus computed in our 
models. This quantity, called is obtained by solving 



the equation: 
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In this way we can directly compare with the analysis of 
Ref. |11|] which explicitly indicates various values of the 
parameter as representative of the EOS tested at 

various energies. In the upper panel of Fig. Owe show the 
parameter computed using the chiral-dilaton model 



4 



and we see that it decreases significantly at large den- 
sities, while its dependence on the temperature is rele- 
vant only at small densities. By comparison, the GM3 
parametrization provides a which is roughly con- 

stant at large densities. 

In Fig. m we show that the incompressibility computed 
using the mixed-phase model remains rather large above 
the lower critical density and it becomes really small only 
approaching the upper critical density. It is extremely 
interesting to notice the strong dependence of the incom- 
pressibility on the temperature. For instance, at T=50 
MeV and p = 3po we still obtain K^^ - 220 MeV, but 
at T=90 MeV and p = 3Apo we obtain I'^O MeV 

and at T=100 MeV and p = 3.6po then R-^ - 120 
MeV. The density-temperature region where the incom- 
pressibility is extremely small is actually rather limited, 
because for T > 150 MeV and p > 4.5po the pure quark 
matter phase is reached, for which the incompressibility 
is large again. Notice also that, while at lower tempera- 
tures the values of the incompressibility is rather large, 
when T exceeds ~ 120 MeV the incompressibility is so 
small that, for practical purposes, the pressure can be 
assumed to be constant. 

A possible scenario based on the mixed-phase forma- 
tion is the following. At 40A GeV the collapse of the 
elliptic flow has been interpreted as a signal of a first 



order phase transition to quark-gluon plasma [1, . At 
lower energies, a mixed phase is produced, but due to 
the existence of two conserved charges the pressure is 
not constant and therefore dramatic manifestations of 
the formation of a mixed-phase are not expected. For 
instance, no wiggle in the behavior of the directed flow 
vs. rapidity is expected, since the appearance of wiggles 
is related to a strong reduction of the incompressibility 
[111. Only rather sophisticated analysis, as e.g. the ex- 
traction of the bulk modulus, can reveal the formation 
of a mixed phase. It will be interesting to test in future 
experiments if the oscillatory behavior of the directed 
flux, observed at 40A GeV, does gradually take place at 
energies exceeding ~ lOA GeV as suggested by our 
analysis. 

In conclusion, a significant softening of the EOS can be 
obtained either via chiral symmetry restoration, if large 
densities are reached, or via the formation of a mixed 
phase of quarks and hadrons. The gradual reduction of 
the bulk modulus discussed in our analysis applies to all 
situations in which the mixed phase forms at not too 
large energy densities, as e.g. in the scenarios discussed 
in Refs. 
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